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ABSTRACT 

We report detection of four IRAC sources in the GOODS-South field with 
an extremely red color of H— [3.6]>4.5. The four sources are not detected in 
the deep HST WFC3 H-band image with H;jmjt=28.3 mag. We find that only 3 
types of SED templates can produce such a red H— [3.6] color: a very dusty SED 
with the Calzetti extinction of Ay=16 mag at z=0.8; a very dusty SED with the 
SMC extinction of Av=8 mag at z=2.0~2.2; and an IGyr SSP with Ay ~0.8 at 
z=5.7. We argue that these sources are unlikely dusty galaxies at z<2.2 based on 
absent strong MIPS 24/im emission. The old stellar population model at z>4.5 
remains a possible solution for the 4 sources. At z>4.5, these sources have stellar 
masses of Log(M,,/M0)=lO.6~11.2. One source, ERS-1, is also a type-II X-ray 
QSO with L2^8A:eV=l-6xlO^'' erg s~^. One of the four sources is an X-ray QSO 
and another one is a HyperLIRG, suggesting a galaxy- merging scenario for the 
formation of these massive galaxies at high redshifts. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies:formation 
— infrared: galaxies 

1. Introduction 

Extremely Red Objects(ERO) are of great interests to modern astrophysics. A simple 
red color criterion generally selects two types of galaxies: those at high redshifts and those 
with heavy dust extinction. With rapid progress in telescope apertures and detectors, this red 
color selection always leads to new types of galaxies or galaxies at record-high redshifts. After 
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large format near-infrared array cameras became availa ble for astronomica l surveys, people 



started to detect galaxies with very red R— K colors (lElston et al.l Il988l EROs, R— K>5 
or I— K>4). EROs were so rare in the early days that they were thought to be abnormal 
objects at very high redshifts. There have been more and more EROs detected by larger 
aperture telescopes with more advanced IR array cameras. Most EROs with R— K>5 are 
now identified as elhptical and dusty galaxies at 0.6 < z < 1.5 ( iThompson et al.l Il999 



McCarthv et al. 


2001; 


Cimatti et al. 


2002) 



spectroscopically i dentified as a dusty Ultra-Luminous InfraRed Galaxy (ULIRG) at z = 1.44 



flElbaz et al.ll2002[ ). This source is analog to a local ULIRG, Arp220. ISmail et al.l (120021 ) 
suggested that most dusty EROs at high redshifts are LlRG/ULIRGs. The Spitzer IRAC 
permits very fast imaging of sky in mid-infrared bands with great depth. IWilson et al.l (|2004[ ) 
found that 17% of their IRAC 3.6/im selected sample are EROs at z>l. 

Red color criteria are practically diversified, and applied to almost all kinds of photome- 
try in optical and IR bands. But the physics for this type of criteria are limited to following: 
(1) the Lyman break at 9 12 A; (2) the Balmer and the accumulated absorption line breaks at 
3648A and 4000A; or (3) dust extinction. Red color caused by the Lyman/Balmer break can 
be used to estimate redshifts. In most deep broad band imagin g surveys, one coul d not tell if 
a red color is due to Lyman/Balmer Break or dust extinction ( ISteidel et al.ll2003l ). An addi- 
tional color in longer wavelength bands is usually applied together with the red color criterion 
to select galaxies at high redshifts. For example, U— g and g— R colors were used to select 
galaxies at z=3 where the Ly man break shifts bet ween U and g bands, commonly known as 
U drop-out for red U— g color (ISteidel et al.ll2003l ). The d rop-out technique w as applied in 
much longer wavelength bands to select galaxies at z=6~ 9. IPranx et al.l J2OO3I ) used the NIR 
color J— K>2.3 to select Distant Red Galaxies (DRGs) with the strong Balmer /4 000A break 



shiftin g in between the J and K at z~2. The NIR spectroscopy for DRGs by iKriek et al. 
( I2OO7I ) shows that about half of their sample are passive evolved galaxies at z~2, and the rest 



are dusty galaxies in a much wider redshift range. Several groups idenitfied 24/im luminous 
and optically faint or invisible sources with R— [24]>14.2 {f2il Ir > 1000) as very dusty 
ULIRGs at z~2. These sou rces are confirmed spectroscopically by Spitzer IRS and ground- 



based optical spectroscopy (IHouck et al.ll2005l : lYan et al 
20091 ). 



20071 : iDev et al.ll2008l : iHuang et al. 



In this paper we report detection of four galaxies with extremely red colors of H— [3.6]>4.5 
in the GOODS-South field. Only one similar source, a submillimeter galaxy (SMG) called 
GOODS 850-5 (aka GNIO) in the GOODS-North field, was ever found to have H-[3.6]>4.5. 
This SMG was also detected by the Submillimeter Array (IWang et al.l 120071 ) with a high 
angula r resolution of ~2" , permitting identification of its counterparts in shorter wavelength 



bands. IWang et al.l ( 120091 ) performed ultra-deep J and H band imaging for this source with 
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NIC3 camera on HST. A total of 16 orbits of HST observation, reaching a nanoJansky depth 
in the F160W band, yields no detection for this source. Based on this red H— [3.6] color, they 
argue that its redshift is at z=4^6.5. La ter, detection of CO (4-3) from this source confirms 
its redshift at z=4.05 ( iDaddi et al.ll2009[ ). This study provides the first look at properties of 
this new type of object. More sources of this kin d will be detected i n the Cosmic Asembly 
Near-infrared Deep Extragalactic Legacy Survey (IGrogin et al.ll201ll . CANDELS). 



2. Deep IR Imaging of GOODS-South 

The Great Observatories Origins Deep Survey(GOODS) is the d eepest mult i -wave length 



survey with space telescopes including HST, Spitzer and Chandra (IDickinson II2004J ). The 



depth of GOODS IRAC 3.6/im imaging reaches sub-microJansky level. Th e deep NIR imag- 



ing of the GOODS-South field was selected for the Early Released Science (lO'Connell 112010 
ERS) for the Wide Field Camera 3 (WFC3), a fourth-generation UVIS/IR imager aboard 
HST. We construct an H-selected multi-wavelength catalog including YJH+IRAC photome- 
try in the ERS covered region. The IR images have very different angular resolutions: 0.03" 
for the HST WFC3 YJH band images and ~2" for the Spitzer IRAC 3.6-8.0 fim images. A 
photometry program called TFIT is specifically designed to perform photometry on a lower 
resolution image with input information of obje ct positions and light distributions measured 



in a high resolution image ( iLaidler et al.ll2007| ). The TFIT program convolves a PSF ker- 



nel to the high angular resolution stamp image for each object to construct lower angular 
resolution image templates and fit them to the lower angular resolution image. In our case, 
we ran the TFIT to perform photometry on the IRAC 3.6/im image for the H-band selected 
galaxies detected in the ERS F160W image. The TFIT also produces a residual image af- 
ter subtracting all H-band detected galaxies in the 3.6/im image. We visually inspected the 
residual image and found four IRAC sources detected at 3.6/im with no H-band counterparts 
shown in Figure [TJ The limiting magnitude for the input H-selected sample is H=28.3 mag 
at 3(7 level, therefore these sources are fainter than H=28.3 mag and have colors redder than 
H-[3.6]=4.5. 

We searched for counterparts of these sources in all available wavelength bands in the 
GOODS-South field. All four sources are detected in the remaining 3 IRAC bands. None of 
these sources is detected in the GOOD-South ACS BVIZ images with the 5cr limiting mag- 
nitudes of 28.65,28. 76,28. 17. and 27.93 respectively. The K-band is the only band available 
in between H and 3.6/im, permitting further constraint on its SED and photon aetric redshift 



The 5a limiting magnitude for the K-band image of the GOOD-South field ( iRetzlaff et al. 



20101) is 24.4 mag and none of our sources is detected in K band. FIR observation is also 
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critical in determining properties of these red sources (IWang et al.l l2009f ) . ERS-3 is clearly 
detected at 24/im and ERS-2 is marginally detected at ~3(j level. We inspected the Herschel 
SPIRE deep imaging of the CDFS, and found only ERS-3 is marginally detected at 250/im 
and 350/im. The PSF for the SPIRE 500 /xr a band image is to o broad ( 36.6") to permit 
accurate extraction of flux density for ERS-3 (IHuang et al.ll201ll ). ERS-3 is also detected at 
1.4gHz with fi.4gj^2=29.2±8/iJy. The remaining three sources are not detected in radio with 
fiAgHz <24/i Jy. No submillirneter/millimemter source is detected in the locations of these 
four sources JScott et al.ll2009yWei^7aDl2009h. Another source, ERS-1, is an X-ray source 
in the Chandra 2Ms catalog (lAlexander et al. 1 l2003h . ERS-2 and ERS-4 are detected only in 
4 IRAC bands (Table.l). 



3. SEDs, Photometric Redshifts, and Properties of the Extremely Red 

Objects 



For three out of the four sources in this study, only NIR-I-IRAC flux densities are 
available for their photometric redshift estimation. The most predominant feature in their 
SEDs is the extremely red color of H— [3.6]>4.5. We first rule o ut that those sourc es are 
brown dwarves. A brown dwarf with T=600K has H— [4.5]>4.0 (ILegget et al.ll2010h . such 
a brown dwarf should have [3.6] — [4.5]>1 due to the methane absorption at 3.6/im in its 
photosphere. All our sources have [3.6] — [4.5] <0. 5. It is unlikely that this red color is due 
to the Lyman break at z>15. Galaxies with a strong Balmer /4000A break at 3<z<8 can 
have very red H— [3.6] colors. Recently iRichard et al.l (120111 ) detected a lensed source at 
z=6.02 behind cluster A383 with H— [3.6] =1.5, arguing that this is a possible passive galaxy. 
A few more lensed galaxies w ith extreme red op tical-MIR color are identified to be dusty 
galaxies at either z~2 or z>6 ( iBoone et al.ll201ll : ?). Their H— [3.6] colors are only in range 
of 0.5<H-[3.6]<2. 



3.1. SED Fitting 



We model the H— [3.6] color using stellar population models of both BC03 ( jBruzual fc Chariot 



20031 . BC03) and the upgraded model CB07 fICharlot fc Bruzualll2007n emphasizing on AGB 



star contribution in the rest-frame NIR bands. The templates are constructed with vari- 
ous stellar populations of the solar metallicity and a very wide range of dust extinction of 
0<Ay <25. The star formation history used in constricting the template set includes single 
burst, exponential decreasing with various e-fold times, and constant rates. Two types of 
extinction curves are used in the SED templates: the widely used Calzetti extinction curve 
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for galaxies (jCalzetti et al.ll2000l ) and the SMC extinction curves (iGordon et al.ll2003l ). Both 
extinction curves are only up to 2.2/im, while our detected photometry points for the four 
sources are in 3.6< A <8.0/im. We ext end both curve to the IRAC bands using the MIR 
dust extinction curve in 3.6< A <24yU,m (jChapman et al.l 120091 ) . 



We first fitted our model templates to GOODS 850-5 to investigate what kind of stellar 
population and how much amount of dust extinction can make such a red H— [3.6] color. 
GOODS 850-5 is already known at z=4.05, thus provides better constraint on stellar pop- 
ulation and dust extinction. Both BC03 and CB07 models with either Calzetti or SMC 
extinction yields a similar result for GOODS 850-5 : IGyr old single stellar population with 
modest extinction of Av'=2.4~3.6. The best fit template is a n IGy r single stellar population 
model with the Calzetti extinction of Ay=3.6. IWang et al.l (120091 ) obtained a similar model 
template for their best fitting but yielding much higher redshifts at z=6.9. 

We argue that the four objects in this study are at the same redshifts: they have very 
similar SEDs, and their positions are very close to each other, with a mean distance of ~1.5' 
to their closest neighbors. We fitted the SED templates to the six IR photometry points 
(H, K, and 4 IRAC bands) for each object in the sample. Our fitting yields two extreme 
solutions with the Calzetti exinction: a very dusty template with Ay=16~18 at z~0.8 and an 
old stellar population template with z~5.7 and Ay ~0.8 (Figure [2]). The templates with the 
SMC extinction yields a similar dusty solution with Av'=7~8 at z~2.2. By applying heavy 
dusty extinction with Ay >7 to templates, its shape and the resulting photometric redshift 
are only determined by extinction curves. For example, the SMC extinction curve yields a 
photometric redshift of Zp=2.2 for our objects, which is caused by a dip at 1.25/im in the SMC 
extinction curve f A ( 1 . 65Mm) / Ay =0 . 1 69 , A(1.25/im)/Ay=0.131, and A(0.81/xm)/Ay=0.567, 
Gordon et al.l (120031 )). With a very high Ay value, this feature is amplified. At z=2.2, this 
extinction dip shifts to the IRAC 3.6;um band to make H— [3.6] redder and [3.6] — [8.0] bluer. 
The photometric redshifts obtained with the SMC extinction curve are mainly driven by this 
feature. 

There are generally three final solutions (Figure [2]) in our SED fitting for the 4 objects: 
a dusty template at z=0.8 with the Calzetti extinction of Ay=16; a dusty template at z=2.2 
with the SMC extinction of Ay=8; and an old stellar population template with age of IGyr 
and Ay <1. Though each solution has a shghtly different minimum x of 3< Xmin <6, we con- 
sider each solution equally possible. All three SED models are able to produce H— [3.6]>4.5, 
and require extreme conditions in the galaxies: either extremely dusty of Ay >7 or even 
Ay=16 or very massive galaxies at z>4.5, both of which are very rare in current extragalactic 
surveys. 
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3.2. Extremely Dusty Galaxies at z<3? 



In the first solution, a galaxy witli tlie Calzetti extinction of Ay ~16 at z=0.8 can liave 
H— [3.6]>4.5. At z=0.8, tlie IRAC 3.6/im band probes tlie rest-frame K-band. Our sources 
have 3.6//m flux density of f3.6=0.6~1.5/iJy, implying that their stellar masses are less than 
5x10^ Mq. Most galaxies with such a small stellar mass at z=0.8 are blue galaxies with no 
dust extinction. M82 is a dusty galaxy with lower stellar mass of 4x10^ M^rj , with heavy 
dust obscuration (5<Av' <51) only occurring in its center ( iBeirao et al. 1 120081 ) . The whole 
M82 appears much bluer than our objects. The H— [3.6] co lor at z=0.8 is equivalent to the 
rest-frame I— K color. M82 has I— K=0.82 (jPale et al.l 120071 ). because most stars in the disk 
of M82 are in the outside of its dusty region. Thus an object like M82 at z=0.8 would be 
detected in the ERS H-band imaging. This solution, however, requires that the whole galaxy 
should be in heavy obscuration. Only ULlRGs have such a obscured morphology. Using M82 
central region SEDs, we predict that the MIPS 24/im flux densities for the first 3 sources 
would have f24=50~70/iJy, well above the FIDEL MIPS 24yum limiting flux density. Only 
ERS3 is marginally detected at 24/im with f24 =39/xJy, the rest are not detected at 24/im. 
We argue that this scenario is least possible. 

The second solution is the template with the SMC extincti on of Ay=7~8 at z=2~2.2. 



Ther e are many Dust Obscu red Galaxies(DO G) identified at z~ 2 ( Houck et al.ll2005 : Dev et al. 



2008 : iBussmann et al.1120091 ) . Several groups (IHouck et al.ll2005l : lYan et al.ll2007l : iHuang et al. 



2OO9I ) performed mid-infrared spectroscopy for DOGs and detected a very strong silicate ab- 
sorpti on feature at 9.8/im in their spectra, indicating a very heavy dust extinction. iBussmann et al, 
( 120091 ) took high resolution H-band imaging of these sources in the Bootes field with NICMOS 
on HST to study their rest-frame optical morpholog ies. We identified these DOGs in the 
Bootes IRAC photometry catalog (lAshby et al.ll2009[ ). and found that DOGs were generally 
red with 1.5<H— [3.6]<3.3, and luminous at 3.6/im with f3.6=5~60/xJy. The four sources 
in this study are much fainter at 3.6Atm, and have a mu ch redder color of H— [3.6]>4.5. 
On the other hand, the DOGs in IBussmann et al.l (120091 ) have a 24-to-8/im flux ratio of 
f24/f8=40~350. If the 4 sources were indeed fainter DOGs at z~2, they should have a 24/im 
flux density of f24 >120/iJy, much higher than the 24/im limiting flux density in GOODS- 
South. Based on much redder H— [3.6] and fainter 24/im flux, we argue that these objects 
are at higher redshifts than the DOGs at z~2. 



3.3. Massive Galaxies at z>4.5? 



The IGyr SSP template with the Calzetti extinction of Ay=0.8 at z~5.7 can also fit 
the SEDs of the four objects, very similar to the best-fit template for GOODS 850-5. In this 
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scenario, the red H— [3.6] is mainly due to the Balmer/4000A jump shifting in between H 
and 3.6/xm bands at z>4. Figure |3] shows that the old stellar template fits to SEDs of the 4 
objects. The resulting photometric redshifts have a large error of 0.4< o"(z)<l.l. We argue 
that the 4 objects are at the same redshifts. Thus by adopting the best a{z) in the 4 objects, 
these objects should be at z>4.5 at 3a level. 

ERS-3 is detec ted at 250 and 350 /xm, thus has a very strong FIR emission, very similar 
to GOODS 850-5 fauang et al.lboilh . The best fit SED models for both GOODS 850-5 
and ERS-3 are old stellar population models. We propose a two-component SED model 
to reconcile the old stellar population and FIR emission in these objects: an old stellar 
population and a very dusty star-forming component. The star- forming component is so 
dusty, sinaliar to those dusty galax ies detected at z~2 ( iHouck et al.l l2005t lYan et al.l 120071 : 
Dey et al.ll2008l : iHuang et al.ll2009l ). that its optical/NIR SED is dominated by the old stellar 
population component. For example, a dusty component with the same stellar mass and 
Av=Q at z>5 would only contribute 10% increase at 8 micron, and much lower percentage 
in the shorter IRAC bands. Assuming a typical dust temperature of Tdu.st=4:0K and redshift 
of z=5.7, we calculate the FIR luminosity for ERS-3 as Log(Lj7/^/LQ)=13.1. The FIR-to- 
Radio flux ratio is a bout q=2.26, consistent with q values for submillimeter galaxies at z > 4 
(IHuang et al.ll201ll ). The remaining 3 objects are not be detected by Herschel SPIRE, but 
may still be IR luminous galaxies with just Log{L fir/I-'q) <13.1. 

The SSP model fitting also yields stellar mass of Log(M*/M0)=lO.6~11.2 for our sources 
(Figure 13]). Spectroscopic confirmed galaxies at z~5.7 including both Lym an-break (LBGs) 



and Ly-a (LAE) galaxies have a typical stellar mass of Log(M^/M (7))<10 (lYan et a. 



Lai et al.l 120071 : lYounger et al.l l2007l : iRichard et al.l l201ll ) . Recently iMarchesini et a. 



2006 



(2010) 



argued that very massive galaxies were already formed at 3<z<4. Theoretically, iLi et al. 
( 120071 ) argued that QSOs at z>6 resided in a massive halo of M~8x10^2Mq, and stellar 
mass for their host galaxies can be as high as IO^^Mq. We have another piece of evidence 
consistent with these syst ems being massive. E RS-1 is also an X-ray source detected in 
the Chandra 2Ms survey (lAlexander et al.ll2003l ). It has x-ray flux densities of fo.5-2keV = 

6.50 ± 1.50 X lO-^^er^ s'^ cm-^. The 
=1.6xl0^^er5f at Zp = 5.7 assuming 



7.55 ± 2.14 X I0~^'^erg s"^ cm"^ and f2-skev = 
hard X-ray Luminosity for this source is lj2-skev 
no absorption correction. Its X-ray-to-optical-fiux ratio {fx/ fR) is higher than 60 and the 
hardness ratio is ~1, thus ERS-1 is an obsc ured type-II QSO. A typical black hole mass 
for su ch a QSO at z=3~6 is 10^ ~ 10^° M© jNetzer 



2003 



Shemmer et al. 2004: Fan et al. 



20061 ). Assuming a typical Ha FWHM of 2000 km/s for a QSO, we convert th e L';i_st.pv to 
black h ole mass for ERS-1 as MBff=5 xlO^ using the relation proposed by ISarria et al. 



mm . iTrakhtenbrot fc Netzerl (120101 ) argued that a host galaxy with 10^ M© black hole 
has a typical stellar mass of M* ~10^^ M0, consistent with the stellar mass we derived for 
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ERS-1. 



4. Summciry and Discussion 

We identified four IRAC sources in tlie GOODS-Soutli field with extremely red color 
of H— [3.6]>4.5. The only known source with a similar H— [3.6] color is GOODS 850-5, a 
SMG in the GOODS-North field. We argue that the four sources must be at the same 
redshift based on the following facts: they have similar rest- frame optical/NIR SEDs; and 
they are spatially very close to each other with a mean angular distance ~ 1.5'. Only 

3 types of templates can produce H— [3.6]>4.5: a very dusty template with the Calzetti 
extinction of Av=16 mag at z=0.8; a very dusty templates with the SMC extinction of 
Av'=8 mag at z=2.0; and an IGyr SSP model with Ay ~0.8 at z=5.7. By comparing 
the 4 objects with local dusty galaxies and DOGs at z~2, we argue that they are unlikely 
dusty galaxies at z=0.8 or z—2.2 based on absent strong 24/xm emission. The old stellar 
population model at z>4.5, with the best fit at z=5.7, remain a possible solution for the 

4 sources. One of our sources, ERS-3, is also detected by Herschel at 250//m and 350//m, 
yielding Log(Li^/ij/L0)=13.2. We propose a two-component SED model for these sources: 
an old SSP component dominating their optical-to-MIR SEDs and a very dusty star-forming 
component mainly contributing to their FIR SEDs. The SED fitting yields stellar masses 
of Log(M^,/MQ) = 10.6~11.2 for the four sources. One source, ERS-1, is also a type-Il x-ray 
QSO with L2-8fceV=l-6x 10'^'^er(7 s~^. Based on the MBH-^buige relation for high-z QSOs, 
ERS-1 should have a massive bulge of Log(M=^/M0)=ll. One of the four sources is an X- 
ray QSO and another one is a HyperLIRG, suggesting a galaxy-merging scenario for the 
formation of these massive galaxies at high redshifts. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407, and with the NASA/ESA HST obtained at the Space Telescope Science 
Institute, which is operated by the association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. 

Facilities: Spitzer(IRAC), HST(STIS), CXO(ASIS). 
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Fig. 1. — Stamp images for the 4 H-band drop-out objects in the GOODS-South field in 
the F160W and 3.6/im bands. ERS-3 is marginally detected at 2 4, 250, 350Atm, and 20 cm. 
ERS-1 is an X-ray source detected in the Chandra 2Ms imaging (lAlexander et al.ll2003l ). 
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Fig. 2. — The likelihood contours as a function of redshift and dust extinction Ky. for the 
best-fit SED of ERS-1. SED fitting for the remaining objects yields the same solutions. The 
left panel is the contour with the IGyr stellar population model and the Calzetti extinction 
and the right panel with the SMC extinction. The IGyr stellar population model with the 
Calzetti extinction of Ay=3.6 is also the best fit for GOODS 850-5 at z=4.05. 
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Fig. 3. — Optical-to-MIR SEDs for the H-band drop-out sources in ERS. The best-fit tem- 
plates to the SEDs are dusty SSP models with E(B-V)=0.2~0.35 and Zp = 5.7. The fitting 
also yields stellar mass of Log(M^,/M0)=lO.6~11.2 for the 4 sources. We also plot two dusty 
model templates against observed SEDs. 
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Table 1. Infrared flux densities for the H-band Dropout sources 



Name 


RA 


DEC 


F160W 


K 


3.6/im 


4.5^m 


5.8/im 




24/im 


ERS-1 


53.084726 


-27.707964 


0.0066±0.002 


0.2904±0.1398 


1.23±0.03 


1.97±0.04 


3.04±0.24 


3.25±0.26 


-11.7±3.8 


ERS-2 


53.132749 


-27.720144 


0.0036±0.002 


-0.0431±0.1681 


1.54±0.20 


1.56±0.10 


2.03±0.26 


3.21±0.28 


14.7±4.0 


ERS-3 


53.060827 


-27.718263 


0.0019±0.002 


0.0727±0.1398 


1.05±0.03 


1.41±0.05 


1.75±0.23 


2.24±0.25 


39.5±4.3 


ERS-4 


53.167161 


-27.715316 


0.0054±0.002 


0.2628±0.1382 


0.57±0.06 


0.66±0.06 


1.25±0.25 


0.87±0.28 


5.9±3.6 



Note. — 



All flux densities in this table are in unit of iJ,3y. 



